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ABSTRACT 

A method for ray-tracing through n-body simulations has been recently pro- 
posed. It is based on a periodic universe covered by simulation boxes. Photons 
move along appropriate directions to avoid periodicity effects. Here, an improved 
version of this method is applied to simulate lensed CMB maps and maps of lens 
deformations. Particle mesh n-body simulations with appropriate boxes and res- 
olutions are used to evolve the nonlinear inhomogeneities until present time. The 
resulting maps are statistically analyzed to look for deviations from Gaussianity. 
These deviations are measured -for the first time- using correlations for config- 
urations of n directions (n < 6). A wide range of angular scales are considered. 
Some interesting prospects are outlined. 

Subject headings: cosmic microwave background — cosmology: theory — large 
scale structure of universe 



1. Introduction 

In this paper, lensed and unlensed maps of the Cosmic Microwave Background (CMB) 
temperature distribution are simulated and analyzed. We are particularly interested in es- 
timating deviations with respect to Gaussianity. Any calculation in this field requires the 
previous assumption of both a background universe and a distribution of cosmological in- 
homogeneities. Taking into account results from the analysis of the WMAP (Wilkinson 
Microwave Anisotropy Probe) first year data (Bennet et al. 2003), the cosmological back- 
ground is assumed to be an inflationary cold dark matter universe with cosmological constant 
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and reionization. The reduced Hubble constant is h — 10~ 2 H = 0.71 (H being the Hub- 
ble constant in units of Km s~ 1 Mpc~ 1 ), the density parameters corresponding to baryons, 
dark matter, and vacuum, are Qf, = 0.04, = 0.23 and Q\ = 0.73, respectively, the total 
density parameter is Q = + + — 1 (inflationary flat universe), and the matter 
density parameter is Q m = Qb + &>d = 0.27. A total reionization is assumed at redshift 
z — 17. Cosmological tensor perturbations are absent, whereas the scalar ones are adia- 
batic energy fluctuations with a Zel'dovich spectrum and a Gaussian distribution. Under all 
these assumptions, the CMBFAST code (Seljak & Zaldarriaga 1996) has been used to get 
the angular power spectrum of the CMB in the absence of lensing (primary anisotropy plus 
the integrated Sachs- Wolfe effect and reionization imprints; hereafter dominant anisotropy) 
and then, a method based on the Fast Fourier Transform (Saez et al. 1996) is used to sim- 
ulate unlensed Gaussian maps of the temperature contrast A D corresponding to dominant 
anisotropy. These maps do not appear to be fully Gaussian because they are generated 
and analyzed by means of discretized numerical processes. Deviations from Gaussianity are 
studied below in an appropriate interval of small angular scales. The simulated maps of 
dominant anisotropy are deformed by lensing. 

The CMB is lensed by cosmological objects evolving in the linear, mildly nonlinear, and 
strongly nonlinear regimes. All these structures deviate the propagation direction of the 
CMB photons. The resulting deviations deform the unlensed maps to produce lensed ones. 
Linear and mildly nonlinear structures can be studied by using analytical or semi-analytical 
methods and, consequently, lens deformations caused by this kind of inhomogeneities can 
be calculated without numerical n-body simulations of structure formation; however, the 
lens effect of strongly nonlinear structures, as galaxy clusters and substructures, is often 
studied by using ray-tracing through n-body simulations. The history of the researches 
about lensing and its status at the end of the last decade can be found in Jain, Seljak, & 
White (2000); in that paper, the authors gave very exhaustive bibliography and described 
and used the classic ray-tracing method based on multiple plane projections. Afterwards, 
another interesting method was proposed by White & Hu (2001). This second method is 
based on many independent simulations, with different boxes and resolutions, tiling the 
photon trajectories. Finally, a third method using an unique simulation has been recently 
proposed (Cerda-Duran et al. 2004). It is used here after some improvements. The main 
features and advantages of this last method are now pointed out. 

The use of a n-body simulation constraint us to work in a periodic universe tiled by 
identical boxes. Photons move in this unrealistic periodic universe and the effect of pe- 
riodicity is important. It magnifies lens deviations. This magnification depends on the 
observation direction and, fortunately, there are directions leading to negligible periodicity 
effects. They are hereafter called preferred directions. In order to understand the existence 
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of these directions (minimizing periodicity effects) let us consider the two-dimensional (2D) 
sketch displayed in Fig. 1, where squares tile a bi-dimensional universe, O is the observer, 
circles are clusters, and lines are photon trajectories. If photons move parallel to the square 
edges (dashed line), they enter neighbouring squares through the same point and pass near 
the same structures in all the squares; hence, lens deviations caused by repeated structures 
add and, consequently, a wrong magnified accumulative lens effect arises. Now, the question 
is: what happens if photons follow the solid line? In such a case, photons enter successive 
boxes through different points separated a distance L as that represented in Fig. 1. For ap- 
propriate values of the angle 0, distance L becomes greater than: (i) the distance at which 
clusters produce relevant deviations of the photon directions (a few Megaparsecs) , and (ii) 
the spatial scales of significant cluster correlations (~ 30 Mpc); hence, distances L > 30 Mpc 
suffice to ensure that photons are deviated by distinct independent clusters located in dif- 
ferent uncorrelated regions. Since there are structures more extended than clusters in the 
n-body simulation, a cutoff avoiding all the spatial scales larger than a scale L cut ~ L is per- 
formed in the peculiar potential given by the n-body simulation. It is worthwhile to notice 
that the simulation uses all the involved scales to calculate the gravitational potential and to 
evolve particles. The scales greater than L cut are subtracted from the resulting gravitational 
potential only to estimate lensing. Using the proposed cutoff, CMB photons are lensed by 
distinct structures in successive boxes and, furthermore, too large scales which are not well 
described in our simulation box are eliminated. Finally, it is worthwhile to remark that, as 
a result of periodicity, there are no discontinuities in the points where photons pass from a 
box to the next one (hereafter called crossing points); as discussed in the next paragraph, 
this fact can be important in order to get a good estimate of lens deviations. 

Roto-translations of the squares could be used to avoid the effects of periodicity. If 
photons move through the dashed line (parallel to the square edges) but a translation of 
length L is performed when photons arrive to the crossing points; then, photons are lensed by 
different clusters in distinct boxes (as it occurs for the preferred directions in our procedure); 
nevertheless, roto-translations produce unavoidable discontinuities in the crossing points. It 
is due to the fact that the photon arrives to an edge at a certain point P, but it enters the 
next square -after translation- through a different point Q. The same occurs for rotations 
and roto-translations and, also, when we use the tiling method proposed by White & Hu 
(2001), in which, different simulations are performed in each square. These discontinuities 
are numerical artifacts and its importance cannot be neglected a priori. 

If the 3D case is considered, the situation is analogous; then, our universe is covered 
by boxes and the trajectories parallel to edges are not appropriate to estimate true lens 
deviations minimizing periodicity effects; nevertheless, if an appropriate observation direction 
is chosen and a cutoff is consistently performed, a large number of boxes can be crossed 



-4- 



through different independent regions without reentering in the zone where the CMB photon 
moved initially. Let us discuss this fact in detail. A Cartesian coordinate system is considered 
in this discussion; its origin O is placed at the center of a box, and the axis are parallel to 
the box edges whose length is L^ ox . The observer is at O and n is the unit vector in 
the direction of a line of sight. This vector can be written either in terms of spherical 
coordinates 9 and or in the form n = (cos a , cos j3 , cos 7) , where a, (3, and 7 are the 
angles formed by n and the x, y, and z-axis, respectively. This second notation is the most 
appropriate for us. Evidently, if a CMB photon go through the x-axis a distance OX, then, 
it advances the distances OY = (cos (3/ cos a)OX and OZ = (cos 7/ cosa)OX along the y 
and z-axis, respectively. Hence, if the photon crosses n boxes in the direction of the x-axis, 
it crosses m = (cos/3/ cosa)n (p = (cos 7/ coso;)n) in the direction y (z). On account of 
these considerations, preferred directions are easily obtained; as an example, for quantities 
n = 5, m = 1, p = 1 + (1/6), and L = 256 Mpc, if the photon crosses a box face orthogonal 
to the x-axis at point M and it crosses the successive parallel face at point N, the distance 
MN in the y-z plane is L = Lb ox [{m/n) 2 + ((p — Vj/n) 2 ] 1 ! 2 = 51.9 Mpc and, furthermore, 
the photon only return to the starting point after crossing n/(p — 1) = 30 boxes in the x 
direction; namely, after going through a comoving distance of ~ 8000 Mpc in the preferred 
direction, which is the distance from z ~ 5.2 to z = in the background universe under 
consideration. At redshift z = 5.2 the box size (256 Mpc) subtends an angle of ~ 1.83 
degrees. Since the observer is located at the center of a box, for a 1.83° x 1.83° map, the 
initial photon positions at z — 5.2 -involved in our ray-tracing procedure- are all located 
inside the same initial box and on a sphere having a comoving radius of ~ 8000 Mpc. For 
greater maps, initial positions cover various boxes and, consequently, the resulting maps 
would have the imprint of the spatial periodicity. That is true even if each photon crosses 
independent regions through a direction close to the preferred one. A large enough number 
of small 1.83° x 1.83° maps allow us a good estimate of correlations for angles smaller than 
~ 0.25° (£ greater than ~ 800). 

Units are chosen in such a way that c = 87rG = 1, where c is the speed of light and G the 
gravitation constant. Whatever quantity "A" may be, A L and A stand for the A values on 
the last scattering surface and at present time, respectively. The scale factor is a(t), where 
t is the cosmological time, and its present value, a , is assumed to be unity, which is always 
possible in flat universes. 
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2. Formalism 

Since lensing produces changes in the propagation direction of the CMB photons, the 
temperature contrast A observed in the n direction is the dominant temperature contrast 
A D corresponding to another direction n ( A(n) = A D (n ) ). The difference 5 = n — n is 
the deviation due to lensing; hence, 

A(n) = A D (n + 6) (1) 
The deviation field due to cosmological structures is (Seljak f996): 

5 = -2 I ° W(\)V ± (p dX , (2) 

where V^0 = — n An A V0 is the transverse gradient of the peculiar gravitational potential, 
and W(X) = (X L - X)/X L . The variable A is 

The integral in the r.h.s. of Eq. (2) is to be evaluated along the background null geodesies. 

The peculiar potential involved in Eq. (2) can be decomposed in the form = 
0i + 02, where 0i (02) is the potential created by the scales smaller (greater) than L cut . For 
appropriate values of the separation scale L cut , potentials 0i and 02 are independent and, 
furthermore, 0i can be obtained with a n-body simulation in a box having a size larger than 
L cut (a few hundreds of Megaparsecs) , whereas the potential 02 can be simulated in a very 
big box (a few thousands of Megaparsecs) by using the potential of the linear approximation. 
There are two kinds of boxes tiling the universe, which are crossed by the same null geodesies. 
Lcut = 60 Mpc is a good choice as it is discussed below. The total deviation given by Eq. 
(2) is the addition of the deviations corresponding to the independent potentials 0i and 2 . 

In the potential approximation (Martinez-Gonzalez et al. 1990), which is valid for linear 
inhomogeneities no larger than the horizon scale and also for nonlinear structures, function 
satisfies the equation: 

A0 = -a 2 (p-p m0 ) , (4) 

where p m0 = fl m oPcrit is the background energy density for matter. In the case of clusters, 
substructures and scales smaller than L cut , the peculiar potential 0i is given by an appro- 
priate n-body simulation -based on Eq. (4)- at each time step. As it is well known, the 
potential 2 created by linear scales larger than L cut has a constant spatial profile and it 
is proportional the ratio Di(a)/a, where -Di(a) is the growing mode of the scalar density 
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fluctuations in the model under consideration (see Peebles (1980) and also Fullana & Saez 
(2000)). At present time, the spatial scales producing significant lens correlations at angular 
scales a < 0.25° (£ > 800) are not larger than the horizon scale (3000/i _1 Mpc) and, conse- 
quently, either Eq. (4) or its Fourier counterpart can be used to get <p2o', namely, to get the 
constant spatial 02 profile. For any of the above potentials, the integral (2) is performed on 
the background null geodesies 



passing by the point O where our observer is located and having directions n. 

The chosen n directions point towards the region of the sky we are mapping and cover 
the map in an appropriate way. We always use a regular coverage. These directions define the 
map nodes where temperature contrasts are calculated and, consequently, they fix the map 
resolution and its size. Once the lines of sight are chosen, unlensed CMB maps of dominant 
anisotropy A D (n) (hereafter U maps) are built up by using the method described by Saez et 
al. (1996) (see §1); hence, the directions n point toward the nodes of the U maps. The lens 
deviations of the potentials 0i and 2 are calculated for the same directions (map nodes) by 
means of Eq. (2). These partial deviations are added to get the total deviation, which is used 
to obtain the direction n + S; finally, the exact Eq. (1) is used -without any approximation- 
to calculate a lensed A(n) map (hereafter L map). Since the direction ft + 5 does not point 
-in general- towards a node of the U map, appropriate interpolation methods are necessary 
to get the quantity A D (n + 5) involved in Eq. (1). The U maps to be lensed has a resolution 
close to one arcminute {£ ~ 10000) and, consequently, interpolations are performed on scales 
much smaller than the angular scales affected by lensing (800 < £ < 3000, see below). The 
high resolution of the U maps allow us the use of linear interpolation (which does not produce 
deviations from Gaussianity) , and also nonlinear interpolation based on splines, which would 
only alter Gaussianity for £ ~ 10000. Finally, the map of the differences A(n) — A D (n) is 
also calculated (hereafter D map); it is the difference (node by node) between a L maps and 
its associated U map; hence, relation L=U+D has an evident meaning. 

Either the n-body potential or the linear one can be used to calculate the integral (2) 
in position space, but we can also work in momentum space; in fact, as it was detailed in 
Cerda-Duran et al. (2004), after simple Fourier algebra, one easily get: 



x % = x l + X(a)ri 



(5) 





where 





with k± = k — (ft ■ k)n and B 
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The n-body simulations are performed using a PM code (Hockney & Eastwood 1988) 
which was used, tested, and described in detail in Quilis et al. (1998). This code gives 
function 5^ at each time step and, then, for a given vector n, the integral of the r.h.s. of Eq. 
(7) can be performed to get function F^(n). Once this function has been calculated, Eq. (6) 
is a Fourier transform which gives the deviation S(x ,n) corresponding to the direction n 
and to the observer located at point x Q ; hence, the deviations associated to all the observers 
located in the nodes of the Fourier box in position space are simultaneously calculated for 
a given direction n. This information is useful to estimate correlations as it was analyzed 
in Cerda-Duran et al. (2004). Similar procedures were used by Aliaga et al. (2002) to 
study the Rees-Sciama effect in momentum space. In the case of linear scales larger than 
L cut , the values of 5^ necessary to calculate the integral in Eq. (7) are given by the linear 
approximation to the evolution of density perturbations. 

Our U, L and D maps must be statistically analyzed. In order to do that, correlations 
associated to n directions (or n points on the last scattering surface) are estimated for n < 6. 
Given a map of the variable ( and m directions, the angular correlations are: 

C m = (C(ni)C(^)---C(n m )> , (8) 

where the average is over many maps. Some of these averages are estimated, in next section, 
for U, L and D maps and for m < 6. Correlations depend on the relative positions of the 
chosen directions, namely, they depend on the figures that these directions draw on the last 
scattering surface. In this paper we have used the sets of directions displayed in Fig. 2, 
where the basic correlation angle, a, is defined as the angle formed by the two directions 
of the case n = 2. In Gaussian statistics it is well known that (Peebles 1980): (i) all the 
C m correlations vanish for odd m and, (ii) for even m, all the correlations can be written 
in terms of C 2 . Hereafter, these even Gaussian correlations are denoted C gm (the suffix g 
stands for Gaussian). For the configurations of Fig. 2, one easily obtains the relation: 

C g4 (a) = 2C\{a) + C%(V2a) , (9) 

leading to C gi (0) = 3Cf (0) for a = 0, and the equation: 

C ge (a) = 3C 2 3 (a) + 2C 2 (a)C 2 2 (v / 2«) + 2C 2 (v / 5a)Q 2 (a) + 
+ AC 2 (a)C 2 (2a)C 2 (V2a) + 2C 2 (v / 5a)C|(v / 2«) + 
+ C 2 (a)C 2 2 (V5a) + 2C 2 (a)C 2 2 (2a) , (10) 

which takes on the form C ge (0) = 15Cf (0) for a = 0. Furthermore, according to Bernardeau 
(1997) and Kesden et al. (2002), if the unlensed maps are assumed to be Gaussian, the non- 
Gaussian lensed ones have vanishing C m correlations for odd m values. Our simulated U 
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maps are almost Gaussian by construction, with small deviations from Gaussianity associated 
to numerical generation; hence, only these numerical deviations could lead to small odd 
correlations without any physical meaning. We are not interested in them; Anyway, the 
correlations C3 and C5 have been estimated for the configurations of Fig. 2 and the resulting 
values have appeared to be compatible with zero. On account of these considerations, the 
simplest method to estimate deviations from Gaussianity is as follows: (a) extract the C2, C4, 
and C 6 correlations from the simulated maps, (b) calculate C g4 , and C g6 functions from C 2 
assuming Gaussianity, namely, using Eqs. (9)-(10) and, (c) compare C 4 and C 6 correlations 
with functions C 9 4 and C 9 q, respectively. 

Now, let us describe the working method used to calculate lens deviations caused by 
the potential 0i (scales smaller than L cut ). First a preferred direction is fixed and a set of 
neighbouring directions are chosen to have a uniform pixelisation of a small map of about 
2° x 2°. This is the crucial point of the method. Afterwards, either Eq. (2) or Eqs. (6)-(7) 
are solved to get deviations and, finally, these deviations are used to deform U maps of 
dominant anisotropy. Now, let us look for appropriate values of the box size and resolutions 
to be used in our PM n-body simulations (0i calculation). The angular power spectrum 
of the lens deformations is expected to be particularly relevant for I > 800; this means 
that correlations must be calculated for a < 0.25° and, consequently, the angular size of 
the simulated maps should be greater than ~ 1°. Our map is chosen to be the image of a 
box face orthogonal to the line of sight of its center and located at the initial redshift; with 
this choice, photon directions do not spread on many boxes at a given time, and there is 
no any undesirable effect produced by the spatial periodicity of the fictitious universe under 
consideration. For this type of maps, the angular size depends on both Lb ox and z in ; it can 
be easily proved that, in our background universe, for Zi n = 5.2 and Lb ox = 256 Mpc, the 
resulting map has a large enough edge of 1.83°. Furthermore, as stated in §1, the separation 
distance L between the regions crossed by CMB photons in successive boxes is 51.9 Mpc, 
the number of boxes crossed between redshifts 5.2 and 0. is ~ 30, and there are hundreds 
of cluster in the box. Three types of simulations with different resolutions are used in this 
paper: (i) Low Resolution (LR) simulations involving 128 cells per edge and 128 3 particles, 
(ii) Intermediate Resolution (IR) simulations using 256 cells per edge and 256 3 particles, 
and (iii) High Resolution (HR) simulations with 512 cells per edge and 512 3 particles. It is 
evident that even our HR simulations, with a cell size of 0.5 Mpc and an effective spatial 
resolution of two cells are not good enough to describe the high density contrasts of core 
clusters; nevertheless, as it will be discussed in detail in next section, higher resolutions 
seem to be unnecessary to study the lens deformations of the CMB maps. Fig. 3 shows 
the column density for particular LR (top), IR (middle), and HR (bottom) simulations. 
Evidently, structures have developed in the three cases, but they have reached different 
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levels of evolution. The diffuse figure of the top panel (LR simulation) indicates a low level 
of evolution (extended structures). The other panels exhibits a better definition because 
structures have evolved more. 

A distance L = 51.9 Mpc suffices to deal with clusters and substructures (see §1), but 
a cutoff eliminating all the scales greater than a maximum one L cut (close to L or smaller 
than it) is necessary to ensure that photons cross different large scale structures in distinct 
boxes. Cutoffs at scales L cut of 60 Mpc, 30 Mpc, 15 Mpc and 7.5 Mpc are independently 
used in next section. A few words about these cutoffs and their consequences are worthwhile. 
The n-body simulations use all the scales corresponding to the chosen values of Lb ox and A 
(spatial resolution); nevertheless, once the n-body simulation has been developed up to a 
certain time, any of the mentioned cutoffs can be performed in the peculiar gravitational 
field responsible for the lens effect (all the scales larger than L cut are thus erased). The 
n-body particles evolve under the forces associated to the full field inside the box, whereas 
the final potential -after the cutoff- is only used to calculate lens deviations using Eq. (2) 
(or Eqs. (6) and (7)). 

In order to perform a certain cutoff, the density contrast, 5g, and its Fourier transform, 
8^ are calculated using n-body simulations and, then the Fourier transform of the total 
peculiar potential inside the box, </>£, is calculated using the Fourier counterpart of Eq. (4): 
^ ^k/k 2 . Afterwards, quantities 0g are avoided for any k < In / L cut and, finally, an inverse 
Fourier transform is used to get the peculiar potential 0i after cutoff. Scales of 60 Mpc, 
30 Mpc are larger than the scale 8/i _1 Mpc ~ 11.3 Mpc and, consequently they are linear, 
but the scales 15 Mpc and 7.5 Mpc are mildly nonlinear. All these scales are well described 
in a box having 256 Mpc per edge. 

Any cutoff should produce two main effects: (i) it should erase too large spatial scales 
which are not well described in the n-body simulation, and (ii) the unaltered scales should 
be small enough to ensure that the CMB photons cross distinct structures in different boxes 
(along the preferred directions). The erased scales are either linear or mildly nonlinear and, 
consequently, they can be studied either using the linear approximation to gravitational in- 
stability (Bardeen 1980) or applying appropriate techniques as the Zel'dovich approximation 
(Zel'dovich 1970), the adhesion model (see Sandarin & Zel'dovich (1989) for a review), the 
lagrangian perturbation approach (Moutarde et al. 1991), and so on; hence, lens deforma- 
tions produced by scales greater than L cut can be calculated without n-body simulations. 
Let us now consider the cutoff at L cut = 60 Mpc. Some scales smaller than 60 Mpc are 
nonlinear and couple among them and with other scales, the n-body simulation is a good de- 
scription of the coupled evolution of all these scales. Once the lens effect produced by scales 
smaller than 60 Mpc has been estimated, the effect produced by strictly linear scales greater 
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than 60 Mpc (complementary scales) can be analytically or numerically calculated using the 
linear approximation and then, taking into account that linear scales do not couple among 
them, and also that these large scales are not expected to be significantly affected by the 
nonlinear ones (which are much smaller), the resulting effect of the complementary scales can 
be added to the effect of scales smaller than 60 Mpc to get the total lens deformations. For 
L C ut values much smaller than 60 Mpc, e.g. 7.5 Mpc, all the complementary scales (greater 
than 7.5 Mpc) are not linear and, furthermore, nonlinear scales smaller and greater than 
7.5 Mpc are coupled; these facts lead to a complicate situation. By these reasons, highly 
linear L cut values, e.g. 60 Mpc, are preferable in order to get the total lens effect, but mildly 
nonlinear values as 7.5 Mpc can be used to get interesting information (see next section). 

In this paper, the lens effect produced by scales smaller than L cut is calculated from 
an initial redshift z in ~ 5.2 to present time. This value of z in seems rather arbitrary, 
nevertheless, in next section it is proved that, for the scales under consideration, the lensing 
produced from decoupling to z — 5.2 can be neglected. 

For L box = 128 Mpc, the size of the resulting maps is a little small (~ 0.9°), the number 
of crossed boxes seems to be too large ( ~ 60), the L value ( ~ 25 Mpc ) is only marginally 
admissible to deal with independent clusters in neighbouring boxes (see §1), and the cutoffs 
must be smaller than ~ 25 Mpc; hence, this box size (and similar ones) seems to be a 
little small. There is only one advantage for L box = 128 Mpc: that the cell size of the 
n-body simulations can be reduced up to 0.25 Mpc, with the same computational cost as in 
the HR simulations -in boxes of 256 Mpc- used in this paper. In the case of large values 
Lbox > 256 Mpc, a cell size as small as 0.25 Mpc requires too much memory; fortunately, 
such a high resolution seems not to be necessary to deal with CMB lensing (see §3). On 
account of these considerations, it seems that L box = 256 Mpc is a good choice to obtain 
appropriate maps allowing statistical analysis. 

Finally, a few words about the calculation of the deviations produced by <p2] namely, by 
the scales greater than L cut = 60 Mpc. In this case, a big box of 5000 Mpc is considered with 
a resolution of ~ 10 Mpc (512 nodes per edge). The potential 2 created by the linear scales 
greater than 60 Mpc is calculated, at present time, by using the Fourier counterpart of Eq. 
(4) and the power spectrum P(k) of the model under consideration (with the normalization 
og = 0.9) and, then, this potential is evolved taking into account that it is proportional 
to Di(a)/a. The lens effect produced by scales larger than L cut = 60 Mpc is calculated 
from decoupling (z in ~ 1100) to present time. Three of these boxes cover all the photon 
trajectories from this initial redshift to present time. For an appropriate preferred direction, 
the distance L is close to 2000 Mpc and, consequently, we can say that, for the scales 
which produce significant contributions to the lens correlations we are estimating (£ > 800), 
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photons cross different boxes through distinct linear inhomogeneities. In this case, there is no 
problem to get a map greater than those of 1.83° x 1.83° associated to the n-body simulations. 
We can easily obtain a map of deviations having a surface sixteen times greater (hereafter 
a 7.32° x 7.32° map), which can be superimposed to a mosaic of sixteen 1.83° x 1.83° maps 
of (j>\ deviations to get a large map of total deviations. 

3. Results 

Two types of maps (hereafter 1 and 2) have been created and analyzed. Type 1 maps 
are small simulations of the partial lens effect produced by scales smaller than L cut ; among 
them we can distinguish: U (unlensed), L (lensed) and D (deformation) maps. All these 
maps have 128 x 128 nodes and 1.83° x 1.83° (excepting two cases described below). Lens 
deviations can be calculated with LR, IR and HR simulations and the associated maps are 
denoted LR-L, LR-D and so on; hence, we have seven sets of maps: the first one is a set of 
eight hundreds U maps to be deformed by lensing, each of the remaining sets (LR-L, LR-D, 
IR-L, IR-D, HR-L and HR-D) contains sixteen thousands maps generated by deforming the 
U maps with the deviations of twenty n-body simulations. Each of these seven sets has 
been analyzed to get both the averaged angular power spectrum and the mean of the C m 
correlations for 2 < m < 6 (see §2). Type 2 maps correspond to the lens effect produced by 
all the scales. Only U (unlensed) and L (lensed) maps (no D maps) are considered in this 
case. All these maps have 512 x 512 nodes and 7.32° x 7.32°. The U maps are lensed by 
total lens deviations, in other words, they are lensed by the superimposition of a map of the 
deviations produced by structures with spatial scales larger than L cut , and a mosaic tiled 
by sixteen 1.83° x 1.83° small maps of lens deviations produced by scales smaller than L cut . 
The discontinuities of the mosaic in the edges of the small maps do not significantly affect 
the estimate of correlations for angular scales a < 0.25°, which are much smaller than the 
size of these small maps. 

L maps are deformations of the U maps and, consequently, neither L and U nor D and U 
maps are statistically independent. According to the relation L=U+D (see §2), the following 
equation is satisfied: 

(L,L)-(U,U) = (D,D) + 2(D,U) . (11) 

The l.h.s of this equation is the difference between the C 2 correlations of the L and U maps. 
The corresponding Cn quantities have a strongly oscillatory character as it was emphasized in 
Seljak (1996) and Hu (2000) (see also Fig. 5). The same methods used to get the oscillatory 
form of the solid line of Fig. 5 (from L and U maps), have been also used to get the Cg 
coefficients of the D maps (associated to the first term of the r.h.s.). The resulting coefficients 
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do not oscillate (see e.g. Fig. 9). The oscillations associated to the l.h.s. of Eq. (11) are due 
to the cross correlations of D and U maps appearing in the second term of the r.h.s. 

We begin with the analysis of small U and D maps of type 1. This study is performed 
with the essential aim of testing the method designed in §2 to look for deviations from 
Gaussianity, this method is used in two oposite situations: almost Gaussian U maps and fully 
non-Gaussian D maps. For U maps, the question is: can we use the method of §2 to detect 
small deviations with respect to Gaussianity due to both the map making method and the 
correlation analysis? The D maps have been considered to verify that large deviations from 
Gaussianity are easily detected and measured with the same method. For U maps, results are 
presented in the left panels of Fig. 4. Right panels show deviations from Gaussianity obtained 
from the HR-D maps. In this Figure, solid lines display correlations directly extracted from 
the simulated maps, whereas dashed lines show quantities associated to the C2 correlation 
function -extracted from the simulated maps- in the case of Gaussian statistics (C S 4 and 
CgQ in Eqs. (9) and (10)). In the top left panel, the resulting C 2 (a) correlation is shown 
in an appropriate a interval where lens correlations are significant (see the top right panel). 
The middle (bottom) panels display the ratio 3Cf (a)/C,4(a) (15Cf (a) / 'Cq(o)) . The values of 
these two ratios are unity for a vanishing correlation angle (a = 0). Small relative deviations 
of the dashed (Gaussian statistics) and solid (map statistics) lines of the left panels suggest 
that the U maps are rather Gaussian by construction. A part of these deviations is associated 
to the use of small maps. Using more extended maps and a greater total coverage, these 
deviations decrease because of a smaller sampling variance (a better estimate of C 2 , C 4 , C 6 , 
C S 4 and C g o). For a large enough coverage, the residual deviations from Gaussianity would 
be essentially due to map making (see below). Right panels are the averaged correlations 
obtained from the HR-D maps. These panels show deviations between solid and dashed lines 
which are much greater than those of the corresponding left panels; hence, D maps are not 
Gaussian at all. The ratios presented in the bottom and middle panels are very useful to 
make well visible -in the figure- the deviations from gaussianity of the U maps (left panels). 
In the right panels, dashed curves present a singular point where the C2 correlation of the 
D maps vanishes (see the top right panel). Very similar results are obtained from LR-D and 
IR-D simulations. Once it has been verified that, even for small type 1 maps, the method of 
§2 works, it can be used in other situations. 

Let us now analyze L and U maps of type 2. The differences between the correlations 
of the L and U maps are hereafter called correlation increments. L maps (not D maps) are 
simplified simulations -without noise, contaminants and so on- of what should be observed 
in experiments and, consequently, correlation increments (excesses or defects with respect 
to the theoretical correlations in the absence of lensing) play an important role in order to 
study whether lensing imprints can be observed in a given experiment. An accurate numer- 
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ical calculation of the difference between the angular power spectra of L and U maps {Ce 
increments) is not trivial. There are various numerical problems. In general, the methods 
designed to extract the Ce coefficients from small maps lead to the worst results in the £ 
intervals where the curvature of the angular power spectrum is great; namely, close to its 
maxima and minima. The extracted Ce coefficients deviate from their true values in these 
intervals and the smaller the map the greater these deviations. Furthermore, Ce increments 
are differences between two similar angular power spectra (which are calculated with numeri- 
cal techniques leading to small errors) and, evidently, such a subtraction -although possible- 
could become problematic; however, the Ce coefficients of the D maps (which measure corre- 
lations in maps of pure lens deformations) are not either oscillatory quantities or differences 
between close spectra and, consequently, they can be accurately obtained with small maps 
and used to do many theoretical and numerical studies about lensing (see below). 

The increments of the Ce quantities are exhibited in Fig. 5. As it is well known, the 
spectrum of the L maps is very similar to that of the U maps and, consequently, Ce increments 
are small; see Seljak (1996) and Hu (2000) for previous comparisons of L and U angular power 
spectra. The dotted line of Fig. 5 shows the increments given by the CMBFAST code for the 
model under consideration. The increments of the solid line are calculated in the following 
form: (i) sixteen 7.32° x 7.32° U maps are simulated and, then, these maps are deformed 
with sixteen 7.32° x 7.32° maps of total lens deviations (see above) to get sixteen L maps, (i) 
the Ce coefficients of all the U and L maps are numerically obtained; the method used to do 
that was described in Arnau et al. (2002) and Burigana & Saez (2003). It is a very efficient 
method to analyze small uniformly pixelised squared maps, and (iii) the averages of the Ce 
quantities extracted from L and U maps are subtracted. Solid and dotted lines have the 
same oscillatory character, but the amplitudes of the maxima and minima are not identical, 
it is mainly due to the small coverage of the analyzed maps. In fact, it has been verified 
that the difference between the amplitudes of the dotted and solid lines increases when our 
study is based on less extended maps. Unfortunately, the use of more extended maps has 
a great computational cost and it is not necessary in this paper, where the comparison of 
the dotted and solid lines of Fig. 5 strongly suggest that our numerical methods work. Let 
us now estimate C m correlations for even m values greater than two. It is done with the 
essential aim of analyzing deviations from Gaussianity due to lensing. 

In order to estimate the C m correlations, the sets of directions shown in Fig. 2 are placed 
at N c random positions on the map and, in this way, the average of Eq. (8) is performed. For 
maps with fixed size and resolution, the resulting correlations only depend on: (1) the number 
Ac, and (2) the accuracy of the interpolation method used to get the temperature contrast in 
the directions that do not point towards map nodes. On account of these facts, U and L maps 
of 7.32° x 7.32°, with a resolution close to one arcminute, have been analyzed for different 



-14- 



values of N c and using very distinct interpolation procedures. In this way, it has been verified 
that, for the two interpolations we have tried (bilinear and splines based on four cells), and 
for N c values close to 5 x 10 5 , the accuracies of the resulting C2, C4 and Cq correlations 
are good enough to allow us -using the method of §2- the study and comparison of the 
deviations from Gaussianity in U and L maps. In order to apply that method, fifty maps of 
type 2 are analyzed as follows: (a) the averaged C2, C4 and C% correlations of the fifty maps 
are calculated using a certain interpolation method and a selected N c value, (b) equation 
(9) is used to get C g4 and, (c) the relative differences AC4/C4 = 2(C 4 — C 94 )/(C 4 + C g4 ) are 
calculated. This analysis is done for U and also for L maps. We begin with the U maps. 
Typical relative differences resulting from the analysis of an unique U map are of the order 
of 10 -2 ; whereas the averaged value corresponding to a few tens of U maps is only of the 
order of 10~ 3 . It is due to the fact that, as the number of maps increases, the estimation of 
the correlations is better, and the quantity AC4/C4 (which vanishes in the Gaussian case) 
decreases to approach its residual value due to deviations from Gausianity associated to map 
making (see above). By comparing the averaged values corresponding to thirty, forty and 
fifty U maps, it has been verified that the AC 4 /C 4 ratio is close to 1CT 3 in all the cases. 
The same values of AC 4 /C 4 are obtained for the two interpolation methods and also for 
N c values greater than 5 x 10 5 . The bottom panel of Fig. 6 describes the same study as 
in the top one, but the involved correlations are now Cq and C 9 q. Results are similar. All 
these results strongly suggest that the AC4/C4 and ACq/Cq averaged values corresponding 
to fifty U maps (solid lines of the top and bottom panels of Fig. 6, respectively) measure 
actual deviations from Gaussianity due to map making. Fifty L maps have been obtained 
from the chosen U maps by using independent realizations of lens deviations. These maps 
have been analyzed in the same way as the U maps. It has been verified that for thirty, 
forty and fifty L maps, the resulting AC4/C4 and AC 6 /C 6 values are always around 20 % 
greater than those of the U maps. This percentage is almost independent on the interpolation 
procedure for any N c value greater than 5 x 10 5 . The dotted lines of the top and botom 
panels of Fig. 6 shows AC 4 /C 4 and ACq/Cq relative differences obtained from fifty L maps. 
Independent lens deviations have deformed the U maps in such a way that the AC 4 /C4 
and ACq/Cq quantities of the resulting L maps sistematically exceed those of the U maps. 
These quantities measure the excess of non-Gaussianity in the L maps beyond what is seen 
in the U maps. All these results strongly suggest that, although deviations from Gaussianity 
produced by lensing are only ~ 20 % of those introduced in the U maps by the map making 
procedure, this ~ 20 % can be pointed out with the method of §2. 

Instead of numerical deviations from Gaussianity associated to map making, observation 
maps involve other deviations due to galactic contaminants, the Sunyaev-Zeldovich effect 
(Cooray 2001), and so on, which are coupled to lens deviations from Gaussianity. Before any 
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eventual separation of the different deviations from Gaussianity appearing in observation 
maps: each contribution and its possible cross correlations with the remaining ones must be 
characterized, including lens contribution we are analyzing here. 

Let us now return to the analysis of small type 1 maps with the essential aim of justifying 
that the initial redshift for lensing calculations (Z in = 5.2) and the resolution of the N- 
body simulations have been appropriately chosen. In order to do that, we study both the 
contribution of the different spatial scales to lens deformations of the CMB and the time 
variation of these contributions. Results are presented in Fig. 7. In order to do such a 
study, scales larger than L cut = 60 Mpc (top panels), L cut = 30 Mpc (middle panels), and 
L cu t = 15 Mpc (bottom panels), have been erased in the peculiar gravitational potential 
and, in each of these three cases, the D maps and their Ci quantities have been calculated 
from z in = 5.2 to redshifts 3.9 (pointed-dashed lines), 2.6 (dashed lines), 1.3 (pointed lines) 
and 0. (continuous lines). Left, central and right panels correspond to an arbitrarily chosen 
LR, IR and HR n-body simulation. In all the panels, the angular power spectrum increases 
with the redshift and it is very small at z — 3.9; hence, for the scales under consideration, 
an initial redshift of 5.2 is a very good choice, and negligible lens deviations are expected at 
z > 5.2. Furthermore, the most important part of the effect is produced between redshifts 
1.3 and (compare solid and pointed lines in all the panels). In any panel (fixed cutoff and 
simulation), the maxima of the curves shift to right as the redshift increases. It is due to the 
fact that the angle subtended by the involved scales depends on the redshift. Since this angle 
(£) is a decreasing (increasing) function of the redshift, the maxima shift to right. For a given 
redshift, the curves of the middle panels are shifted to right (left) with respect to those of 
the top (bottom) panels. It is due to the fact that -as a result of the cutoff- the scales of the 
middle panels are smaller (greater) than those of the top (bottom) one and, consequently, 
they subtend a smaller (greater) angle for a fixed redshift. From the comparison of the 
continuous lines in the top and middle panels it follows that spatial scales between 60 and 30 
Mpc, whose effects are (are not) included in the top (middle) panel, are responsible for the 
most important part of the total effect (solid line of the top panel). Similarly, the comparison 
of the solid lines of the middle and bottom panels indicates that the effect produced by the 
scales between 30 and 15 Mpc is greater than that produced by the scales smaller than 15 
Mpc in all the cases. This last effect due to scales smaller than 15 Mpc is very small for 
LR simulations (left bottom panel), but a comparison of the bottom panels -in which the 
n-body resolution increases from left to right- among them lead to the conclusion that this 
effect of scales smaller than 15 Mpc grows as the n-body resolution increases. It is due to 
the fact that the greater the n-body resolution the greater the growing of the scales smaller 
than 15 Mpc. In the case of the HR simulation of the right panels of Fig. 7, a cutoff at 7.5 
Mpc has been also performed. Results are presented in Fig. 8 (same structure as Fig. 7), 
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where we see that the amplitude of the effect produced by scales smaller than 7.5 Mpc is 
much smaller than that produced by the complementary scales (up to 60 Mpc). Of course, 
this effect should increase as resolution does, nevertheless, for the angular scales of interest, 
this effect would be negligible -for any resolution- as compared with that of the remaining 
scales. From the above considerations, it follows that the most significant scales are greater 
than 15 Mpc in all cases. Fortunately, these scales are rather well described by our LR 
and IR simulations and, consequently, no important differences appear among the averaged 
spectra of our LR-D, IR-D and HR-D maps. These spectra are presented in Fig. 9, where we 
see that, for great enough £ values, the larger the n-body resolution the greater the resulting 
Ci quantities. These I values correspond to angles subtended by small enough spatial scales 
which grow more as resolution increases. We see that LR simulations suffice to get a very 
good approach to the angular power spectrum of lens deformations. 

Twenty LR simulations have been used in our study. Each of these simulations leads to 
eight hundred D maps which can be analyzed to get an averaged angular power spectrum; 
in this way, twenty different power spectra appear whose mean has been represented in the 
solid line of Fig. 9. We are now interested in the distribution of these twenty spectra around 
the mean spectrum of Fig. 9. How much separated are these spectra from their mean? In 
order to answer this question, the mean and variance of the twenty Cg, coefficients have 
been calculated for every t value and, afterwards, mean, mean plus variance, and mean 
minus variance have been used to built up solid, dotted and dashed lines of the top panel 
of Fig. 10, which gives a good idea about the dispersion of the spectra associated to each of 
the twenty LR simulations. The middle (bottom) panel of Fig. 10 has the same structure as 
the top one, but it corresponds to IR (HR) simulations. 

Since a LR simulation suffices to get a rather good description of CMB lensing, we have 
developed a new LR simulation in a box of 512 Mpc (256 3 cells and 256 3 particles) using: 
(i) a cutoff at 60 Mpc, (ii) fifteen boxes to cover the photon trajectories from z in = 5.2 and, 
(iii) a preferred direction leading to a separation distance L = 103.8 Mpc. Evidently, the 
situation is better than in previous simulations because of the smaller value of the number 
of boxes and the larger values of L and Lb ox . Since the box size is twice that of previous 
simulations, the new L and D maps are greater than those we have already analyzed, they 
have 3.66° degrees per edge. New U maps have been built up in concordance with the new 
coverage. In the top panel of Fig. 11, we show the spectra of the resulting D maps at redshifts 
3.9, 2.6, 1.3, and 0. (same format as in Figs. 7 and 8). The spectrum is very similar to those 
obtained with previous LR simulations (compare with the top left panel of Fig. 7 taking into 
account that we are considering the spectra of two particular LR simulations). These results 
strongly suggest that boxes of 256 Mpc are large enough to study CMB lensing. 
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Results obtained from Eq. (2), in position space, are now compared with those obtained, 
in momentum space, from Eqs. (6) and (7) (plus statistics based on 128 3 observers, see §2). 
In order to do this comparison, LR simulations as those used along the paper (excepting 
previous paragraph) are considered. The cutoff is performed at 60 Mpc. The main difference 
with respect to previous cases is the size and resolution of the simulated D maps. We have 
not been able to work with 128 x 128 directions -in momentum space- as a result of the 
high computational cost of the method based on Eqs. (6) and (7), which works with 128 3 
observers at the same time. By this reason and with the essential aim of comparing results in 
position and momentum spaces, regular maps covered by 16 x 16 directions with a resolution 
of 5' and an edge of 1.25° have been simulated. Size and resolution are smaller than in 
previous cases and, consequently, only some Ce quantities can be extracted from the new D 
maps with admissible accuracy. In any case, the D maps can be simulated using the two 
methods and identical U maps; afterwards, the Ce coefficients can be extracted using the 
same codes and, finally, results can be compared. The solid line of the bottom panel of 
Fig. 11 shows the spectrum obtained, in position space, by averaging the spectra of twenty 
D maps with the new size and resolution. Dashed and dotted lines exhibit the Q quantities 
found, in momentum space, using only one LR simulation in each case, but considering 
128 3 observers in the statistical analysis (see Cerda-Duran et al. (2004) and §2). Figure 11 
shows that both methods lead to similar curves. If we average over twenty simulations in 
momentum space, the final spectrum is close to that displayed in the solid line of Fig. 11; 
namely, both methods lead to similar final spectra separately. 



4. Discussion and conclusions 

A method for ray-tracing through n-body PM simulations (Cerda-Duran et al. 2004) is 
used to simulate lensed maps of the CMB. Our attention is focused on the study of devia- 
tions from Gaussianity in these maps. Some correlation funct ions as C4, Cg, CgA and Cg% are 
estimated and compared -for the first time- to measure deviations from Gaussianity. Con- 
siderable emphasis is put on the comparison of the lens deformations produced by different 
small spatial scales, and also in the study of the evolution of these deformations. 

For a given map, the relative difference between the correlation functions C4 and C g 4, 
and also between C% and C 9 q measure deviations from Gaussianity. These relative differences 
-AC4/C4 and ACq/Cq- have been calculated for U and L maps. According to our expecta- 
tions, the functions AC4/C4 and ACq/C^ of the U maps take on small values pointing out 
small deviations from Gaussianity due to construction (see §3). This success in the analysis 
of U maps strongly suggest that the same method could be applicable to the analysis of other 
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deviations from Gaussianity (lensing, Sunyaev-Zelovich and so on). In the case of L maps, 
quantities AC4/C4 takes on values which are greater than the values of the U maps. The 
same occurs for AC 6 /C 6 . The differences between the AC4/C4 (AC 6 /C 6 ) quantities of the 
L and U maps are about the twenty per cent of the AC4/C4 (ACq/Cq) values corresponding 
to the U maps. The same values of these differences are obtained for different interpolation 
procedures and whatever the value of the parameter N c > 5 x 10 5 may be; hence, these 
differences measure actual deviations from Gaussianity produced by lensing (see discussion 
in §3). 

L and U maps lead to very similar Ce quantities. Differences between these quantities 
are calculated with CMBFAST and also from the analysis of our simulated maps. These 
differences are displayed in Fig. (5). Both methods lead to comparable results. This fact 
strongly support the ray-tracing procedure, the numerical simulations, and the correlation 
analysis used in this paper. 

Box sizes of ~ 256 Mpc seem to be appropriate for ray-tracing through n-body sim- 
ulations. Much greater (smaller) sizes lead to problems with resolution (L values, number 
of boxes and so on). In these boxes, scales between 15 and 60 Mpc -small enough as com- 
pared with chosen box size- produce the largest part of the lens effect on the CMB, whereas 
scales smaller than 15 Mpc only produces a small part. By this reason, low resolution n- 
body simulations suffice to get a good approximation to the angular correlations produced 
by lensing. For scales smaller than 60 Mpc, the most important part of the lens effect is 
produced between redshifts 2.6 and 0., whereas only a small part of this effect is produced 
at z > 3.9. 

Two independent codes based on distinct approaches have been used to deal with the 
same problem: the estimate of lens correlations produced by scales smaller than L cut . One 
of these codes works in momentum space and involves many observers (Cerda-Duran et al. 
2004) for statistical analysis. In the second code, calculations are performed in position 
space -method applied along this paper- and the observer is unique. Both codes lead to 
very similar angular power spectra after averages. 

Let us now list some prospects. It would be worthwhile: (1) the use of other n-body 
codes (allowing more resolution) with the essential aim of getting more accurate correla- 
tions at small angular scales, (2) the application of appropriate modern methods (wavelets, 
Minkowski functionals and so on) to study deviations from Gausianity in the L and U maps 
we have simulated, (3) the characterization of all the deviations from Gaussianity which 
appear coupled in observation maps, e.g. deviations due to the Sunyaev-Zel'dovich effect, to 
the radiation from the galaxy in the CMB frequencies and so on. The necessary statistical 
analysis could be developed by using either the method of this paper or other methods (see 
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above), and (5) the study of the possible detection of deviations from Gaussianity in maps 
of the PLANCK mission. 
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Fig. 1. — Sketch of a 2D universe covered by squared patches. Circles are structures, the 
observer is located at point O and lines are photon trajectories 






Fig. 2. — Configurations of n directions for 2 < n < 6. The sets of n directions draw these 
Figures on the Last Scattering Surface. 




Fig. 3. — integrated column density along an edge of the simulation box. Top, middle and 
bottom panels correspond to a LR, IR and HR simulation, respectively. 



- 23 - 




0.05 0.10 0.15 0.20 

«[degrees] 




1.0 - 



0.05 0.10 0.15 0.20 

a[degrees] 




0.05 0.10 0.15 0.20 

a[degrees] 



6 



5 



4 



7 







1 1 1 1 1 1 1 1 1 1 1 1 


1 




\ 






r / ^ 






I 1 1 






/ \ 






1 \ 






1 \ 






1 \ 




/ ; 


Z ' \ 
I ' \ 




/ ; 


/ 




/ j 


/ 


\ 

\ / 




/ 


\ / 
\ ' 
\ / 










- = 



0.C5 0.10 0.15 0.20 

a[degrees] 



Fig. 4. — Left (right) panels correspond to U (D) maps. Top, middle and bottom panels 
exhibit functions C 2 , 3C|(a)/C 4 (a) and 15C|(q;)/C 6 (q;) against the correlation angle a in 
degrees. In the solid lines, correlations C2, C4, and Cq are directly extracted from the maps, 
whereas the ratios shown in the dotted lines are calculated using the Gaussian correlation 
C g4: and C g Q defined in the text. The separation between continuous and dotted lines measure 
deviations from Gaussianity. 
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Fig. 5. — Difference ACe in fiK 2 , multiplied by the factor 1(1 + 1)/2tt, v.s. \og£. Difference 
is calculated between the Ce coefficients of L and U maps. Dotted line has been obtained 
with CMBFAST and solid line from the numerical analysis of sixteen 7.32° x 7.32° L and U 
maps (see the text). 
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Fig. 6. — Top panel shows relative differences AC4/C4 between the C 4 correlations extracted 
from the U maps and the C 9 4 Gaussian correlations (see text), v.s. the correlation angle 
in degrees. Bottom panel has the same structure as the top one, but it exhibits relative 
differences AC 6 /C 6 . 
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Fig. 7. — Quantity £(£ + l)Q/27r in /ii^ 2 v.s. log£. Left, central and right panels involve 
LR, IR and HR simulations. Top, middle and bottom panels correspond to L cut values of 60, 
30, and 15 Mpc. Solid, pointed, dashed, and pointed-dashed lines display the C t coefficients 
of the D maps produced by lensing between z in = 5.2 and the redshifts 3.9, 2.6, 1.3, and 0., 
respectively. 
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Fig. 9. — Quantity 1(1+ l)Ce/2n in fiK 2 v.s. \og£. The Cg coefficients involved in the solid, 
dotted and dashed lines correspond to the averaged power spectra of the sets of maps LR-D, 
IR-D and HR-D defined in the text 
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Fig. 10. — Quantity £(£ + l)Q/27r in fj,K 2 v.s. \og£. Top, middle and bottom panels show 
results from LR, IR and HR simulations. Solid lines are the averaged spectra displayed in 
Fig. 9, whereas the dotted and dashed lines show la curves in the spectra distribution (see 
the text for more details). 
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Fig. 11. — Top panel has the same structure as panels of Fig 7. It corresponds to a LR 
simulaton in a big box of 512 Mpc, with a cutoff at L cut = 60 Mpc. Bottom panel represents 
the same function as in the top one, but the Ce coefficients have been obtained from D maps 
with a size of 1.25° and a resolution of 5'. The Ce quantities of the solid line correspond 
to a map obtained from Eq. (2) in position space, whereas in the dotted and dashed lines, 
we present the Ce coefficients of a map, which has been found from Eqs. (6) and (7) in 
momentum space (see the text for more explanations). 



